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ABSTRACT 

Context. The temperature is a central parameter affecting the chemical and physical properties of dense cores of interstellar clouds 
and their potential evolution towards star formation. The chemistry and the dust properties are temperature dependent and, therefore, 
interpretation of any observation requires the knowledge of the temperature and its variations. Direct measurement of the gas kinetic 
temperature is possible with molecular line spectroscopy, the ammonia molecule, NH3, being the most commonly used tracer. 
Aims. We want to determine the accuracy of the temperature estimates derived from ammonia spectra. The normal interpretation of 
NH3 observations assumes that all the hyperfine line components are tracing the same volume of gas. However, in the case of strong 
temperature gradients they may be sensitive to different layers and this could cause errors in the optical depth and gas temperature 
estimates. 

Methods. We examine a series of spherically symmetric cloud models, 1.0 and 0.5 Mq Bonnor-Ebert spheres, with different radial 
temperature profiles. We calculate synthetic NH3 spectra and compare the derived column densities and temperatures to the actual 
values in the models. 

Results. For high signal-to-noise observations, the estimated gas kinetic temperatures are within ~0.3 K of the real mass averaged 
temperature and the column densities are correct to within ~10%. When the S/N ratio of the (2,2) spectrum decreases below 10, the 
temperature errors are of the order of 1 K but without a significant bias. Only when the density of the models is increased by a factor 
of a few, the results begin to show significant bias because of the saturation of the (1,1) main group. 

Conclusions. The ammonia spectra are found to be a reliable tracer of the real mass averaged gas temperature. Because the radial 
temperature profiles of the cores are not well constrained, the central temperature could still be different from this value. If the cores 
are optically very thick, there are no longer guarantees of the accuracy of the estimates. 

Key words. ISM: clouds - ISM: molecules - Radio lines: ISM - Stars: formation - radiative transfer 



> 

(N 

O 

cn 

o 

(N 



X 



1. Introduction 

The dense molecular cloud cores are of particular interest be- 
cause their collapse can lead to the formation of new stars. To 
understand the star formation process, one needs to understand 
the physical and chemical properties of the cores. The tempera- 
ture of the cores, especially as a function of position, is difficult 
to measure and yet it affects everything from core stability to 
chemistry and dust evolution. 

Dust and gas become therma lly coupled only at high densi- 
ties, n(H2)~ 10 cm and above (Goldsmith 2001). Because the 
continuum surface brightness is strongly affected by the warm 
outer cloud layers, observations of the dust emission are not ideal 
to determine the gas temperature at the centre of the dense cores. 
The widths of molecular lines can sometime s be used to derive 
stringent limits on the gas temperature (e.g.. lHariu et al.l [2008) 
but the value of this method is limited by the presence of turbu- 
lent velocity component. 

In observations of most molecules, it is difficult to separate 
the effects of density, optical depth, and temperature. The situ- 
ation is different for ammonia, NH3, whose rotational energy is 
characterised by the quantum numbers (J,K). Because the dipole 
transitions between the different K ladders are forbidden, their 
relative populations depend only on collisions and thus measure 
the gas kinetic temperature. Each rotational state (J,K) is further 
divided to inversion doublets. The inversion transitions (1,1) and 
(2,2) are at 23.7 GHz and 23.72 GHz respectively, and are eas- 
ily observable from the ground. The inversion transitions also 



have a hyperfine structure. For example the (1,1) spectrum con- 
tains 18 components from which at least five separate groups of 
lines are usually discernible. The intensity ratio of the hyperfine 
components provides a convenient tool to estimate the line opti- 
cal depths. Because of these properties, ammonia is a common 
tool for gas te mperature measureme nts, especially in the case 
of dense cores (lHo&Townesl ll983). The high critical density 
(~ 2 • 10 3 cm~ 3 for the (1,1) transition as calculated from the ra- 
tio of the spontaneous and collisional de-excitation rates at the 
temperature of 15 K) guarantees that one is measuring temper- 
ature of dense gas. On the other hand, NH3 is resilient against 
depletion and the abundan ce is observed to drop o nly in the very 
densest and coldest cores (Ber gin & Tafallall2007l) . 

The analysis of NH3 is based on the assumption that the 
{J,K) transitions and their different hyperfine components are 
all tracing the same physical region. Because of the different op- 
tical depths, the weaker hyperfine components could be more 
sensitive to deeper layers in the cloud. The effect should become 
noticeable if the main lines are moderately optically thick and 
the source has strong temperature gradients. 

The temperature structure of quiescent cores can, in princi- 
ple, be calculated from th e balance of the various heating an d 
cooling mechanisms (e.g. lGoldsmithll200ll: TGalli et alJ 12002). 
The gas within the cores is heated by the cosmic rays and the 
external stellar radiation field and it is cooled down by the 
line radiation from a number of atomic and molecular species. 
Depending on the temperature difference between the gas and 
the dust, the gas-dust collisions can result in either heating or 
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cooling of the gas. lJuvela& Ysardl (T20 111) discussed the range of 
possible temperature profiles in dense cores, emphasising the un- 
certainties associated with many of the contributing factors (e.g., 
the external radiation field, the gas-dust coupling, and the abun- 
dances of the species responsib le for the main cooling lines). In 
this paper we use the models of lJuvela & Ysardl d201 lb to exam- 
ine the reliability of the ammonia temperature estimates in the 
presence of radial temperature variations. 

The content of the paper is the following. In Sect. [2] we 
present the methods concerning the core models, the calcula- 
tion of synthetic ammonia spectra, and the derivation of temper- 
ature and column density estimates. The results are presented in 
Sect. [3] where we compare the derived estimates to the corre- 
sponding actual values in the models, in particular studying the 
effects of observational noise and increased column densities. 
The final conclusions are presented in Sect. [5] 



2. Methods 

2. 1 . Core models 



Table 1. The main core models. 



We use the cloud models discussed by Juvela & Ysardl (1201 ll) . 
These are 0.5 Mq and 1.0 Mq Bonnor-Ebert spheres (£=6.5, 
2\in=10K) for which the gas temperatures were calculated un- 
der different hypotheses. This resulted in a set of radial temper- 
ature profiles with varying shapes, the temperature being of the 
order of 10K. While the dust temperature decreases inwards in 
starless cores, the situation can be different for the gas tempera- 
ture because of the efficient molecular line cooling in the outer 
part. In the models, the gas temperature profiles were mostly in- 
creasing towards the centre. When the photoelectric heating was 
included, the situation is reversed and the surface temperature is 
raised close to 20 K. 

The calculations follow Goldsmith (2001) apart from the fact 
that the cooling rates due to line emission are obtained with a 
Monte Carlo radiative transfer program. The cooling rates are 
calculated directly for 12 CO, 13 CO, C 18 0, C, CS, and o-H 2 0. 
The rates of the last two species are multiplied by 10.0 and 
2.0, respectively, to take into account the contribution of other 
species. We examine the following cases without the photo- 
electric heating: the default model ( l def), the low abundance 
model (x), the weak gas-dust coupling model (A g d), the small 
velocity dispersion model (cr v ), and the large velocity gradient 
model (v). The names refer to the temperature calculations in 
Uuvela& Ysardl (l20Tll) . 

The def-model was obtained assuming cr v =1.0kms _1 line 
widths (Doppler width) and constant abundances in the calcu- 
lation of the cooling lines, no large scale velocity field, and a 
coupling betwee n the gas and dus t temperatures that followed 
that presented in Goldsmith (2001). In the cr v -model, the tur- 
bulent linewidth changes as a step function from l.Okms 1 in 
the outer part, outside the half radius, to O.lkms 1 within the 
inner core. The v-model is the only one with a large scale veloc- 
ity field. It has radial infall velocities that increase linearly from 
Okms -1 at the surface to l.Okms -1 in the centre. The name of 
the low abundance model (x) refers to the calculations of the gas 
kinetic temperature where the abundance of the cooling species 
was assumed to be lower at the centre of the core. This was 
implemented as a step function with a factor of ten difference 
between the inner and outer parts. However, in the modelling 
of the NH3 spectra, the fractional abundanc e of para-NH3 is al- 
ways kept at a co nstant value of 10~ 8 (e.g. iHotzel et al. 2001; 
Tafalla et al. 2004). This applies also to the ^-model. The pho- 
toelectric heating is included only in one model (model PEH) 
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where the model cloud is assumed t o be illuminated by an unat- 
tenuated interstellar radiation field (M athis et alj 1983[). For the 
details of the calculations of the 7\i„ profiles, see lJu vela & Ysardl 
(120111) . 



In addition to the above described models, we consider cases 
where the kinetic temperatures are the same but the turbulent line 
width is decreased by a factor of ten in the NH3 calculations. 
In particular, in the central part of the cr v -mod el the linewidth 
then b ecomes practically thermal. We use the lJuvela & Ysardl 
(1201 ll) results only as examples of the possible shapes of the 
temperature profiles and the purpose of the present paper is to 
probe their effect on the analysis of ammonia spectra at a general 
level. 



2.2. Synthetic spectra 

For each model, synthetic para-NH3 (1,1) and (2,2) spectra were 
obtained with the help of radiative transfer calculations. The 
non-LTE radiative tran sfer problem is solved with a Monte Carlo 
program (|Juvelalll997l) . The the model clouds are divided into 
100 shells. To get better resolution in the outer regions where the 
excitation conditions can change more rapidly, the thickness of 
the shells was decreased towards the surface. The radius of the 
innermost sphere and the thickness of the outermost shell are, 
respectively, ~6% and 0.7% of the radius of the model cloud. 

The ammonia spectra are calcul ated following the prescrip- 
tion presented in iKeto et all (120041) . The non-LTE populations 
of the (J,K) states of para-NH3 are obtained from the radiative 
transfer calculations where, within each state, an LTE distribu- 
tion is assumed between the hyperfine levels (i.e., distribution 
according to the statistical weights of the transitions). In the 
radiative transfer calculations, we take into account the hyper- 
fine components of the (1,1) and (2,2) lines. These consist of 
18 and 21 components, respectively. In the final synthetic (1,1) 
and (2,2) spectra, many of these components are blended to- 
gether. The molecula r parameters for the (J,K) lines are from 
the Lamda database ( Schoier et al] 120051) and the frequencies 
and relative int ensities of the hyperfine lines were taken from 
lKukolichl(ll967h . 

The synthetic spectra are calculated for a pencil beam (i.e., 
for a single line-of-sight through the centre of the core) and av- 
eraged over a Gaussian beam with the FWHM equal to the core 
radius. Because the excitation of the molecules varies along the 
line of sight and the hyperfine components have different optical 
depths, the observed intensity ratios between hyperfine compo- 
nents can vary depending on the optical depth and temperature 
structure of the model. Figure[TJshows examples of the simulated 
spectra. 
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Fig. 1. Examples of simulated ammonia spectra. Each frame 
shows the (1,1) (the upper line) and the (2,2) spectral profiles. 
For the plotting, the latter have been scaled by a factor of 3 and 
shifted by -2 K. The text in the upper right hand corner indi- 
cates the model in question and, in parentheses, the turbulent 
line width used in the calculation of the NH3 spectra. 



2.3. Analysis of the spectra 

To derive estimates of the column densities and gas kinetic 
temperatures, we subject the synthetic NH3(1, 1) and NH^(2, 2) 
spectra to standard analysis as descr ibed in iHo & Townes 
( 1983), IWalmslev & Ungerechtsl (1 19831) . and lUngerechts et alJ 
( 1986). The hyperfine structure of the NH3(1, 1) line is fitted us- 
ing a x 2 routine assuming a gaussian velocity distribution and 
that the different hyperfine states are populated according to 
LTE. The fit results are used to calculate the optical thickness 
profile t(v). The peak value r pea k, and the brightness tempera- 
ture in the corresponding channel (usually coincident with the 
intensity maximum rB iPea k) are substituted in the antenna equa- 
tion to derive the excitation temperature, r ex , of the transition. 
The integrated optical thickness, J r(v)dv, where v is the radial 
velocity, and T ex give the column density in the upper transition 
level, N u (l, 1), from the formula 



8tt 1 

NJl,l) = -3- — F(T U , 



^ex) r(v)dv, 



where A u i is the Einstein coefficient for spontaneous emission, 
An = c/vn is the wavelength of the transition, Tn = hvn/ks, 
and the function F(Tq, T) is defined by 



F(T , T) = 



1 



,T /T _ 1 



Adding together the upper and lower levels of the inversion tran- 
sition, the total column density of ammonia in the (1,1) level is 



N(l,l) = N u (l+e 



Tu/T„ 



)■ 



The column density of the (2,2) level, N(2, 2), can be derived 
using two slightly different methods: a) by estimating t and T ex 
from the hyperfine fit as in the case of (1,1); or b) by assuming 
that r ex is the same as for the (1,1) transition and that the (2,2) 
line is optically thin. In the case b the integral J r(v)dv is esti- 
mated from the integrated intensity of the (2,2) main group using 
the antenna equation. Usually only the method b is feasible for 
observational data of cold clouds because the (2,2) satellites are 



very weak and buried in the noise. In the present paper we use 
the method b to do the analysis. 

The rotational temperature, T12, is defined with the 
Boltzmann equation 



2V(2,2) 
N(l, 1) 



-4i.5/r l2 



The kinetic temperature, 7\i n , is estimated using the three 
level approximation (inc luding the levels (J, A")=(l,l), (2,1) 
and (2,2)) as outlined by IWalmslev & Ungerechtsl (119831) . The 
temperature dependence of the ratio of the collisional coeffi- 
cie nts C?.?,-»?.i /C??.^i 1 was estimated adopting the coefficients 
of iDanbv et al.l (119881) . The method is the same as used in 
lHariu et all (Il993l) . The resulting relationship betwe en T12 and 
?\in is very close to that given in the Appendix B of Tafalla et al. 
d2004l) . In cold gas, the sum of the column densities of the (1,1) 
and (2,2) levels is almost equal to the total column density of 
para-NH3 as the populations of the higher rotational levels are 
negligible. We calculate these column density estimates from the 
synthetic observations and compare with the corresponding true 
values in the model clouds. 



3. Results 

3.1. The basic models 

The derived para-NH3 column densities and the gas kinetic tem- 
peratures were compared to the corresponding true values that 
were calculated as the mass averaged quantities in the model 
clouds, averaged over the same beam as what was used in the 
calculation of the synthetic ammonia spectra. 

The main results are presented in Figs. [2]- [4] The profiles 
of the kinetic temperature 7\j n are shown as solid black curves. 
These are a critical input parameter for the simulation of the am- 
monia spectra. Without the photoelectric heating, 7\m increases 
towards the cloud centre. The central temperature is further in- 
creased when, within one half of the cloud radius, the abundance 
of the cooling species is decreased (model x) or the turbulent 
linewidth is decreased. In model A gl d, the gas-dust coupling 
is weakened in the inner core, leading to a small decrease in 
the temperature because 7d us t > Tgas- We c '° not cia i m that all 
these temperature profiles are likely to be found in starless cores. 
However, for the present case the most important fact is that they 
probe a wide range of T^ n profiles that potentially could result 
in a bias in the 7\j n estimates. 

Figure [2] shows the 7\j n estimates for the 1 Mq cloud. The 
horizontal thin and thick lines indicate the true mass averaged 
temperatures for a single line-of-sight through the centre of the 
model cloud (thin line) or weighted with a gaussian beam with 
a FWHM equal to cloud radius. The dashed lines are corre- 
spondingly the 7\j n estimates derived from the NH3 spectra ob- 
served with a pencil beam or the larger gaussian beam. In the 
frames also are given the true and estimated column densities. 
For comparison, Figs. [3] and [4] show the results for 0.5 Mq 
models and for 1 Mq models with a smaller turbulent linewidth 
(o- ^O.lkms- 1 ). 

The temperature estimates derived from the simulated NH3 
spectra are correct to within 0.5 K and the errors of the column 



1 In Juvela & Ysardl j201 1[) the gas-dust cou pling A g d was calculated 
using the formulas giv en in [Goldsmith (2001). The true heat exchange 
rate may be higher (see lYoung et alj2004l) . This would not have a large 
effect on the temperature in the core centres where the gas and the dust 
are well coupled. However, it would modify somewhat the shape of the 
Tkin profiles. 
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Fig. 2. The radial temperature profiles (solid black curves), the 
estimated temperatures (dashed red lines), and the mass aver- 
aged true temperatures (black horizontal lines) for the one solar 
mass models. The thin lines correspond to the data on a single 
line-of-sight through the cloud centre and the thick lines to data 
averaged with a gaussian beam. In the right hand part of each 
frame, the difference between the estimated and the real mass 
averaged 7\in has been scaled by a factor of ten for better vis- 
ibility. The column densities derived from the observations are 
indicated in each frame (units of 10 14 cm -2 ), the subscript C de- 
noting the case with beam convolution. The true column densi- 
ties are given in the first frame. For the v-model, the dotted line 
and the lower text entry indicate the result from a fit with two 
velocity components (pencil beam only). 



density values are ~10% or less. This applies to the observa- 
tions conducted with a pencil beam and most observations with 
a gaussian beam. When one includes the convolution with the 
large beam, the temperature values reflect more the values found 
in the outer part of the cloud. However, compared to the beam 
averaged real kinetic temperature, the values appear to be biased 
towards the T^\ n in the cloud centre. The error is largest for the 
<x v model where it can approach 1 K (Fig. [2j. 

The spectra of the v-model are double peaked because of 
the infall motion that at the centre leads to maximum line-of- 
sight velocities of +lkms~' (see Fig. Q}. The redshifted peak 
is stronger than the blueshifted peak. This is opposite to the 
expected behaviour in contracting clouds and is caused by the 
fact that the temperature is decreasing towards the core cen- 
tre within the innermost ~20% of the cloud radius. When the 
spectra are modelled as the sum of two components, the indi- 
vidual features described with gaussians, the analysis recovers 
the correct solution with good accuracy. In the case of narrow 
lines (<x <0.1 kms -1 ) a single component fit to the two velocity 
components would result in an error of ~ 1 K. On the other hand, 
when that model is observed with a larger beam, only a single ve- 
locity component remains visible. In Fig.|4]the beam convolved 
data were fitted with a single gaussian and the pencil beam ob- 
servations with two components. In both cases the temperature 
estimates remain within ~0.2 K of the correct value. 
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Fig. 3. As Fig. [2]but for the half solar mass core. 
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Fig. 4. As Fig.[2]but with the synthetic ammonia spectra mod- 
elled with a smaller, 0.1 kms 1 turbulent linewidth. For the 
model v, the pencil beam observations are fitted with two veloc- 
ity components and the convolved spectra with a single velocity 
component. 



3.2. The influence of noise 

The influence of observational noise was examined by analysing 
a series of spectra with progressively higher signal-to-noise ra- 
tios. The S/N ratio is defined as the ratio between the peak 
intensity of the (2,2) spectrum and the rms noise per chan- 
nel. The rms noise added to the (1,1) spectra was the same in 
brightness temperature units. In addition to the default model 
(marked 'def'), the models with the temperature distributions 
corresponding to the 'cry' and 'PEH' cases were examined, i.e., 
models with strong inward and outward temperature gradients. 
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Fig. 5. The gas kinetic temperature (plus symbols) and the 
column density (crosses) estimated from spectra with different 
signal-to-noise ratios. For each S/N ratio, the figure shows re- 
sults from five realisations and their average as the dashed lines. 
The horizontal lines indicate the true values in the models. The 
three frames correspond to the default model (def), the cry- 
model, and the PEH model, respectively, of the 1.0 M sun cloud. 



Both the M = 0.5 Mq and M = 1.0 Afo cloud models were 







used, calculating the ammonia spectra with either l.Okms or 
0.1 kms linewidths. Again, the narrower line width is not con- 
sistent with the origin al derivation of the temperature profiles 
(I Juvela & Ysardl 1201 ll) but gives some insight to clouds with 
higher line optical depths. 

Figure [5] shows the estimated kinetic temperatures and col- 
umn densities in relation to their true values in the model clouds. 
The spectra were calculated using a pencil beam. The S/N covers 
a range from ~3 to 300 in multiplicative steps of 1.5. For each 
S/N ratio we show the results from five realizations of spectra, 
the dashed lines corresponding to their mean value. The results 
from the M - 0.5 Mq and the <x v =0. 1 km s _1 were qualitatively 
similar. 



3.3. Models of high optical depth 

In the default 1 M sun model the peak optical depths are T(i,i)=0.7 
and T(2 j 2)=0.013. In the 0.5 M sun cloud these increase to 1.2 and 
0.04, respectively. When the turbulent linewidth was decreased 
from l.Okms -1 to O.lkms -1 , the values were naturally higher, 
4.3 and 0.33 for the (1,1) and (2,2) lines, respectively. However, 
as the satellite lines continue to be mostly optically thin, this is 
not a big problem for the analysis of the spectra. 

To determine when and how the errors increase to signifi- 
cant levels, we took the 1 .0 M sun clouds and calculated a series of 
models where the densities were multiplied with constant factors 
k n between 1 and 64. The temperature profiles are taken from 
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Fig. 6. The gas kinetic temperature (plus symbols) and the col- 
umn density (crosses) when the densities of the 1 .0 M slm model 
are scaled with factors k n in the range of 1-64. The column den- 
sities have been normalized with k n . The solid lines indicate the 
true values and the dashed lines are the values deduced from the 
synthetic observations. The frames correspond to different tem- 
perature distributions as marked in the upper right hand corner 
of each frame. 

the original cr v and PEH models (Fig. |2]i and are not changecQ 
Therefore, the changes will be only because of the increase of 
the density and the column density. Figure [6] shows the result- 
ing 7\i n and N(para-NHi) estimates resulting from observations 
with a pencil beam towards the cloud centre. As expected, the 
estimates show some bias when the column densities are much 
higher than for the original 1 Mq Bonnor-Ebert spheres. In the 
PEH model the 7\i n estimates continue to be very accurate and 
even in column density the error rises to ~ 50% only when the 
density is some 50 times the original values. In this model, the 
real temperature is quite constant within the central regions and 
apparently the optical depths are still not high enough so that the 
spectra would react to the warm but low density surface layers. 
The changes were more dramatic in the cr v -model where the re- 
sults were essentially wrong after k n ~4. In this model the optical 
depths were higher to start with, because of the narrow linewidth 
at the centre. 

Figure [7] shows the actual spectra calculated for k n -\ and 
£ n =32. In the cr v -model, the main group becomes strongly self- 
absorbed and, because the kinetic temperature is higher in the 
centre than at the surface, the intensity of the satellite lines is 
higher than for the main group. This is, of course, a situation that 
cannot be accounted for in the analysis assuming a homogeneous 
medium. 

4. Discussion 



Tafal la etaD d2004l) examined the validity of the kinetic tem- 
perature estimates from NH3 using isothermal models. In the 
present study, we calculate the mass averaged kinetic temper- 
atures in model clouds with realistic temperature profiles. The 



2 We are investigating the effects of different temperature profiles in 
general. However, we note that for high k n the gas kinetic temperature 
should approach the dust temper ature and Tun is be expected to decrease 
towards the core centre (see e.g. Crapsi et al. 2007). The situation could 
change only if the core is heated internally. 
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Fig. 7. Comparison of the original spectra from the one solar 
mass <x v and the PEH models and the spectra obtained after the 
densities have been scaled by a factor of k n =32. For the plotting, 
the (2,2) spectra have been scaled by a factor of 3 (upper frames) 
or 0.5 (lower frames) and shifted by -2 K. 

values are compared to the kinetic temperature estimates de- 
rived with a standard LTE method, using the ammonia spectra 
'observed' from these models. 

The estimated kinetic temperatures were generally found 
to be very accurate. In the case of model clouds with masses 
0.5 Mq and 1 Mq, the difference between the temperature de- 
rived from the NH3 spectra and the real mass-averaged temper- 
ature was of the order of 0.1 K, only rarely approaching 1 K. 
Similarly, the estimated column densities of para-NH3 remained 
correct to within -10%. 

In one of the models the infall velocity produced two veloc- 
ity components whose presence, in the case of large turbulent 
linewidths cry = 1.0 km -1 , was not immediately obvious from 
the line profiles. The analysis completed using a single gaussian 
component still led to only ~ 1 K error. A fit with two velocity 
components resulted in twice as high errors in both the tempera- 
ture and the estimated column density (Fig. [2j». When the turbu- 
lent line width was reduced in the cloud core so that the velocity 
components became clearly visible in the (1,1) profiles, the er- 
ror of a single component fit increased to ~ 1 .5 K and the column 
density was overestimated by 20% (Fig. HJ. However, this time 
the two component fit did recover the correct results with less 
than 10% errors. This increases the general confidence in the re- 
sults of obtained from ammonia spectra observed with adequate 
resolution. It appears unlikely that the effect of hidden features, 
e.g. of multiple velocity components, would be amplified in the 
7\j n and N estimates. 

When the signal-to-noise of the (2,2) spectra is decreased 
below S/N~10, the temperature errors increased in our tests to 
~1 K. The results were not significantly biased so that even for a 
set of sources observed at low S/N, the mean properties should 
be relatively trustworthy. The sensitivity to noise naturally de- 
pends on the width of the lines. 

When the column density is increased much beyond that 
found in our basic models, the (1,1) becomes optically thick and 
there are no longer guarantees of the accuracy of the results. In 
particular, the lines might no longer probe the conditions at the 
centre of very dense cores. We carried out one test with an ad hoc 
scaling of the densities. The behaviour was found to depend very 
much on the line widths (affecting the initial line opacity) and the 
temperature structure of the cloud. In the case with a small tur- 



bulent linewidth, <x v = O.lkrrT 1 s , completely wrong results 
were already encountered with densities a few times higher that 
those of an 1 M sun Bonnor-Ebert sphere. If the fractional abun- 
dance were larger than the 10~ 8 assumed in our models, obser- 
vations of normal hydrostatic cores could be affected. However, 
in starless cores the values are not likely to be much higher 
(iHotzel et al.ll200 lb iTafalla et al.ll200arCrapsi et alJl2007h . 

iTafalla et alJ (120041) found that the excitation temperatures 
of the (1,1) and (2,2) lines varied by a factor of several. This 
variation depends on several factors, the cloud optical depth, 
the volume density, and the variation of the kinetic tempera- 
ture. Figure|8]shows the situation in those of our I.OMq models 
where the difference in the radial kinetic temperature profiles is 
the largest. In the cry model the kinetic temperature increases 
inwards and, in combination of higher volume density and en- 
hanced photon trapping, the excitation temperatures at the cloud 
centre are three times as high as at the cloud boundary. Towards 
the centre also the difference between the excitation tempera- 
tures of the (1,1) and (2,2) transitions becomes more noticeable. 
On the other hand, in the PEH case the outward increasing ki- 
netic temperature helps to keep the excitation conditions almost 
constant. Therefore, it is not surprising that the kinetic tempera- 
ture estimated from the observed spectra remained accurate even 
when the density was increased (Fig. [6]). 

Ammonia was found to be a very reliable tracer of the real 
mass averaged gas temperature within the modelled Bonnor- 
Ebert spheres. The temperature at the centre of a cloud could 
still differ from this number by several degrees. In our models 
this difference was at most ~1 K. However, if the signal from 
very centre becomes weaker because of a lower abundance com- 
bined with a low kinetic and excitation temperatures, the differ- 
ence can be more significant. 

In this paper we have assumed the ammonia abun dance to be 
const ant as a function of the cloud radius (see e.g. ITafalla et alJ 
2006). If the abundance varies, NH3 is likely to provide accurate 
estimates of the average temperature weighted not only by the 
mass but also by the fractional abundanc e. If the ammonia abun- 
dance increases towards the centre (e.g. ITafalla et al.ll2002h the 
result will be a better estimate of the central core temperature. In 
cores denser than considered here the situation may be reversed 
as even NH3 becom es depleted in the very centre of the core 
(lAikawa et al ] l2005h . As an illustration of the associated uncer- 
tainties, we considered the cr v model with 1 .0 Mq, varying the 
abundance between 10~ 75 and 10~ 8 5 as a linear function of the 
logarithm of the radius. One obtains an estimate of T^i n - 10.49 K 
when the abundance increases inwards and 10.59 K when the 
abundance increases outwards. Both are within 0.25 K of the 
value from the previous constant abundance model. Thus, in this 
case the effect of the abundance gradients is not very significant. 
However, we will return to this questions in future papers that 
discuss self-consistent models of the chemistry and of the ther- 
mal balance. 



5. Conclusions 

We have calculated ammonia emission spectra using non-LTE 
radiative transfer modelling and cloud core temperature profiles 
taken from the literature. We have carried out analysis of the sim- 
ulated spectra and compared the results with the true parameters 
of the model clouds. This has led to the following conclusions: 

- For the ~1 Mq hydrostatic clouds, the kinetic temperatures 
and column densities derived from the ammonia observa- 
tions are accurate to better than 10% accuracy. 
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Fig. 8. Comparison of temperatures in the 1 .0 Mq cloud models 
0" v and peh. Shown are the radial profiles of the kinetic temper- 
ature (solid lines) and the excitation temperatures of the (1,1) 
(dashed line) and (2,2) (dash-dotted line). 

- The presence of several velocity components can be identi- 
fied from the line profiles before they have a significant effect 
on the accuracy of the analysis. 

- The uncertainty of the kinetic temperature estimates in- 
creases to ~ 1 K when the signal-to-noise of the (2,2) spectra 
is decreased to S/N~10. 

- In models of higher density, the behaviour depends the ki- 
netic temperature profile of the cloud. Quite wrong results 
can already be encountered with models a few times more 
opaque than the 1 M sun Bonnor-Ebert spheres. 

- Apart from the most opaque cores, the ammonia was found 
to be a very reliable tracer of the mass averaged kinetic tem- 
perature. The temperature at the centre of dense cores could 
still differ from this value by several degrees. 
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